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ABSTRACT 



Context. There are a number of faint compact infrared excess sources in the central stellar cluster of the Milky Way. Their nature and 
origin is unclear. In addition to several isolated objects of this kind there is a small but dense cluster of comoving sources (IRS13N) 
located ~3" west of SgrA* just 0.5" north of the bright IRS13E cluster of Wolf-Rayet and 0-type stars. Based on the analysis of their 
color and brightness, there are two main possibilities: (1) they may be dust embedded stars older than few Myr, or (2) very young, 
dusty stars with ages less than 1 Myr. 

Aims. We present a first K^-band identification and proper motions of the IRS13N members, the high velocity dusty S-cluster object 
(DSO, also referred to as G2), and other infrared excess sources in the central field. Goal is to constrain the nature of these source. 
Methods. The L'-band (3.8yL/m) observations have been carried out using the NACO adaptive optics system at the ESO VLT. Proper 
motions have been obtained by linear fitting of the stellar positions extracted by StarFinder as a function of time, weighted by 
positional uncertainties, and by Gaussian fitting from high-pass filtered and deconvolved images. We also present results of near- 
infrared (NIR) H- and K,-band ESO-SINFONI integral field spectroscopy of the Galactic Center cluster ISR13N. 
Results. We show that within the uncertainties, the proper motions of the IRS13N sources in K^- and L'-band are identical. This 
indicates that the bright L'-band IRS13N sources are indeed dust enshrouded stars rather than core-less dust clouds. The proper 
motions also show that the IRS13N sources are not strongly gravitationally bound to each other. Combined with their NIR colors this 
implies that they have been formed recently. For the DSO we obtain proper motions and a K^-L' -color. 

Conclusions. Most of the compact L'-band excess emission sources have a compact H- or K^-band counterpart and therefore are 
likely stars with dust shells or disks. Our new results and orbital analysis from our previous work favor the hypothesis that the infrared 
excess IRS13N members and other dusty sources close to SgrA* are young dusty stars and that star formation at the Galactic Center 
(GC) is a continuously ongoing process. For the DSO the color information indicates that it may be a dust cloud or a dust embedded 
star. 
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^ 1 ■ Introduction 2006). Currently the "in-situ" scenario appears to be favored by 
r a number of authors (Nayakshin & Sunyaev 2005, Nayakshin 

^ The formation of young and massive stars is an important pro- et al. 2006, Paumard et al. 2006). Also, results by Stolte et al. 

^ cess at the center of the Milky Way (see e. g. Paumard et al. (2007, 2010) seem to rule out the possibility that known com- 

5r 2006, Ghez et al. 2005). In the central half parsec of the Milky pact clusters close to the GC (like e.g. the Arches cluster) could 

Way these objects are organized in at least one disk-like struc- migrate inwards and fuel the young stellar population at the very 

ture of clockwise rotating stars (CWS; Genzel et al. 2003, Levin center. 
& Beloborodov 2003, Paumard et al. 2006) and the existence of 

a second, less populated disk of counter-clockwise rotating stars A prime object to test these scenarios both observationally 

(CCWS) has been proposed by Paumard et al. (2006). It is un- and theoretically is the IRS 13 group of sources. IRS13E (lo- 

clear how these young stars have been formed in the strong tidal cated ~ 3" west and ~ 1.5" south of SgrA*) is the densest stel- 

field of the super-massive black hole (SMBH) at the position of lar association after the immediate vicinity of SgrA* and con- 

SgrA*. There are two prominent scenarios: One includes star tains several massive Wolf-Rayet and 0-type stars (Maillard et 

formation "in-situ" (in an accretion disk; Levin & Beloborodov al. 2004, Moultaka et al. 2005, Paumard et al. 2006, Fritz et al. 

2003, Nayakshin et al. 2006), the other requires the in-spiral of 2010). It is generally considered to be associated with the mini- 

a massive stellar cluster that formed at a safe distance of 5-30 spiral (Moultaka et al. 2005, Paumard et al. 2004). For IRS13E 

pc from the Galactic Center (Gerhard et al. 2001, McMillan & four out of seven identified stars show highly correlated veloc- 

Portegies Zwart 2003, Kim et al. 2004, Portegies Zwart et al. ities (Maillard et al. 2004, Schodel et al 2005) indicating that 
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is probably bound. It is unclear what the origin of the IRS13E 
group of stars is. It is conceivable that such a cluster could have 
been formed in an accretion disk as described by Milosavljevic 
& Loeb (2004) and Nayakshin et al. (2005). However, numerical 
simulations have shown that it is difficult to explain the forma- 
tion of such a dominant feature in a star forming disk (Nayakshin 
et al. 2007). In order to support the cluster in-fall scenario, an 
intermediate-mass black hole (IMBH) was proposed to reside in 
the center of the cluster (Maillard et al. 2004). The existence of 
an IMBH results in a higher dynamical friction and a higher sta- 
bility of the in-falling cluster such that this process in general is 
more efficient. It would allow compact clusters to reach the cen- 
tral parsec of the Galaxy (with >10^Mo) within their lifetimes 
(Hansen & Milosavljevic 2003, Berukoff et al. 2006, Portegies 
Zwart et al. 2006, but see Kim et al. 2004 for a characterization 
of the problems with this hypothesis). However, the presence of 
the IMBH that would be required to stabilize IRS13E is highly 
disputed. Schodel et al. (2005) analyze the velocity dispersion of 
cluster stars. The authors found that the mass of such an object 
should be >7000Mo. However, both the X-ray (Baganoff et al. 
2003) and radio (Zhao & Goss 1998) source at the position of 
IRS13E can be explained by colliding winds of high-mass los- 
ing stars (Coker et al. 2002, Zhao & Goss 1998). These findings 
make the presence of an unusual and massive and dark object in 
IRS13E unnecessary. 

In addition to the IRS13E association there are a number of 
faint infrared excess sources that may be essential for the discus- 
sion of star formation in the central stellar cluster (see Eckart et 
al. 2006a and Perger et al. 2008). In some cases there is spec- 
troscopic evidence for them being associated with young stars 
(Perger et al. 2008). Here we present new proper motion and 
new spectroscopy data on these objects Q. These comprise first 
K^-band proper motion measurements of sources that include the 
IRS13N association, and several sources at a projected distance 
of only a few arcseconds from SgrA*. 

Approximately 0.5" north of IRS13E, a small cluster of un- 
usually red compact sources has been reported (IRS13N; Eckart 
et al. 2004). Muzic et al. (2008) give a detailed analysis of this 
cluster and show that an orbital analysis supports that the clus- 
ter is a comoving group of young stars. A strong infrared excess 
is due to the emission of warm dust (T ~ 1000 K; Moultaka et 
al. 2005). The authors propose two possible explanations for the 
nature of IRS13N: (1) objects older than few Myr and similar to 
bow- shock sources reported by Tanner et al. (2005) or (2) very 
young stars (0.1-1 Myr old). The latter scenario implies more 
recent star formation than what has been assumed so far for the 
GC environment. 

Recently, Gillessen et al. (2012) reported a fast moving 
strong infrared excess source which they interpret as a core-less 
gas and dust cloud approaching SgrA*. Due to its location and 
due to the fact that it is apparently on an elliptical orbit around 
SgrA* we will refer to it in the following as dusty S-cluster ob- 
ject (DSO). The community has started to call the newly found 
fast moving object G2 (e.g. Burkert et al. 2012), reminiscent of 
the gas cloud Gl found by Clenet et al. (2005). As we show 
in the present paper it cannot be excluded that the newly found 
object is rather a star with a gas/dust shell or disk around it, 
than a core-less gas cloud as claimed by Gillessen et al. (2012). 
We therefore propose to call it DSO rather than G2 since the 
acronym DSO (Dusty S-cluster Object) reflects the true nature 
of the source in a more realistic way. While essential proper- 



^ Based on observations collected at the European Southern 
Observatory, Chile 



ties of this source have already been described by Gillessen et 
al. (2012), no identification of it shortward of 3.5/dm has been 
published to date. In other cases of near-infrared excess sources 
(e.g. see Fig. 14 in Eckart et al. 2006a) a detailed investigation 
of their nature had not yet been carried out. We also present first 
K^-band proper motion studies of the DSO. 

After the introduction we describe the observations and data 
reduction in section O In the results section [3] we present the 
derivation the proper motions and then outline the physical prop- 
erties of the sources in sectionlH In section [5] we give a summary 
and the Appendix contains additional maps and graphs. 

2. Observations and Data Reduction 

2.1. NACO near-infrared imaging 

The imaging data in the H-, K^- and -bands were obtained 
with the NAOS/CONICA adaptive optics (AO) assisted im- 
ager/spectrometer (Lenzen et al. 1998, Rousset et al. 1998, 
Brandner 2002) at the UT4 (Yepun) at the ESO Very Large 
Telescope (VLT). IRS7, a supergiant with K^^^ 6.5 - 7.0 mag 
and located -5.5" north from SgrA*, was used to lock the AO 
loop using the NIR wavefront sensor. Observing dates and pixel 
scales for the sources in the central few arcseconds, the IRS13N 
sources and the DSO are given in Tab. [T] Data reduction was 
standard, with sky subtraction, bad pixel correction, and flat 
fielding. Randomly dithered exposures were median averaged to 
obtain final mosaics. The sky background in the H- and K^-band 
was measured on a dark cloud about 400" north and 713" east of 
the GC. In the L'-band, the sky background was extracted from 
a median of the dithered science exposures. These are images in 
which the profiles of single stars are sharply peaked and have 
no broad wings. In the K^-band one can see a diff'raction limited 
stellar image (in some cases with an indication of a diff'raction 
ring) on top of a more extended foot that stems from the resid- 
ual of the AO seeing correction. From these images we selected 
data sets in which the complete NIR PSF close to SgrA* has a 
FWHM of less that 0.25" (estimated through the second order 
moment; optical guide star seeing typically <0.5") and an esti- 
mated FWHM of the diffraction limited part of the PSF of better 
than 1.5 times the diffraction limit (by visual inspection of the 
PSF profiles). The diff'raction limited angular resolutions in the 
Ks- and L'-band are about 42 mas, 56 mas and 98 mas, 
respectively (NaCo instrument manual). Infrared PSF width es- 
timates for all NACO H-band data until 2008 are given with the 
light curves by Bremer et al. (2011). Active optics guide star 
PSF width estimates obtained in the optical for all NACO K^- 
band data until 2010 are given in Fig.20 by Witzel et al. (2012). 
For the years 2002 and 2009-2012 we did not find high quality 
^-band imaging VLT data. 

The sources in the wider field of the central star cluster are 
shown in Fig. [TJ The IRS13E and IRS13N clusters are shown 
in Figs. [2] and O The relative and absolute coordinates of the 
lower left (south-eastern) corner of the shown images are listed 
in Tab. [21 In the IRS13N region 8 sources have been resolved 
and labeled a through k as shown in Figs. IH [51 and [71 The 
nomenclature is taken from Eckart et al. (2004). 

2.2. SINFONI imaging spectroscopy 

Spectroscopy of the region around IRS 1 3 was carried out in June 
2009 using SINFONI, the AO-assisted integral field spectrome- 
ter mounted at Yepun, Unit Telescope 4 of the ESO VLT in Chile 
(Eisenhauer et al. 2003). The AO guiding was carried out in the 
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Fig. 1. Proper motion vectors of dusty objects in the central stellar cluster of the Milky Way (see Tab.©. We show the May 2003 Adaptive 
Optics L'-band image of the central 4" x 4". The arrows indicate the uncertainties in amount (red or yellow portion of the otherwise green arrow) 
and angle of the measured proper motion derived from data between 2002 and 2008 in the H, and L' -bands (see Tab.[T]). 



optical using a R~13.9 star located ~19" from SgrA*. We used 
the K^-band grating providing spectral resolution of 4000; the 
spatial setting of 0.1" per image slice results in a field-of-view 
of 3"x3". The mean Strehl ratio of the cube (measured from 
Bry) is about 7%. The sky was sampled on the same dark cloud 
as in the case of NACO observations, and the total on-source 
time was 14x300 s. Because of the non-quadratic shape of the 
SINFONI pixels (50 mas x 100 mas), the observations were split 
in two blocks observed at position angles of and 90 degrees. 
The two cubes were combined into a single data cube as a final 



step in the reduction procedure. To combine the two data cubes, 
we first de-rotated the second cube by 90 degrees to match the 
orientation of the first one. Both data cubes were re-binned to 
a unique pixel scale of 50 mas x 50 mas using the IDL task 
'congrid' and median-collapsed. The resulting images were then 
cross-correlated to determine the off'sets that can be supplied to 
the pipeline recipe ' sinfo-UtLcube -combine ' . The recipe shifts 
the two cubes and co-adds them to produce the final data cube. 

The reduction and reconstruction of the 3D cubes were car- 
ried out using the ESO supported SINFONI reduction software 
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Table 1. Observational data of the GC used for proper motions 
measurements in H- and K^- and L'-band for the sources given 
in Fig. [T] except the dusty S-cluster object (DSO) for which the 
epochs are given in the text (section [2Jl) . 



package. Bad pixel, cosmic rays, and flat field corrections were 
applied to the 2D raw frames. The 3D cubes were reconstructed 
using calibration frames for the slit-let distances and the light 
dispersion. Intermediate standard-star observations of a nearby 
G2V star were used to correct for strong atmospheric (telluric) 
absorptions. The standard-star observations were done near in 
time and airmass to the target exposures. The G2V spectral 
characteristics were removed by dividing the standard-star spec- 
trum by the well known, high signal-to-noise solar spectrum 
(NSO/Kitt Peak spectra, as provided at the ISAAC/VLT web 
page). 

The flux calibration was done using the known fluxes of the 
IRS13N cluster sources (Eckart et al. 2004; Muzic et al. 2008), 
since we are primarily interested in measuring line ratios. The 
total flux within the SINEONI EOV was scaled according to 
these values. The positional registration of the data cube was 
done via the know positions of the bright IRS13E sources and 
its uncertainty is of the order of 0.05". 



Erom the final reconstructed 3D cube we extracted spec- 
tra and emission line maps. In Eig. [8] and [9] we show the K- 
(intermediate continuum band IB 2.15yum) and L'-band image 
(3.8yum) images with the Bry and [Eelll] -emission lines and the 
CO absorption line maps in contours over-plotted. The contours 




2006 I 1 arcsec = 39 mpc 



Fig. 2. High-pass filtered image of the IRS13E/IRS13N field in K,- 
band. The labels correspond to those in Tab. [3] Source 1 to 6 are sources 
El to E6 in IRS13E. The angular resolution is close to the difiraction 
limit of 60mas. 




Fig. 3. A wavelength integrated 2yt/m SINEONI data cube. Yellow cir- 
cles mark the regions in which the spectra were obtained. The angular 
resolution is of the order of 0.11 arcseconds. 

in Eigs. m [9] and [TOl are labeled with percentages of the peak 
flux. The emission line maps presented in Eig. [8] were created by 
summing up the flux of the emission line and then subtracting 
the average flux density of the continuum to the left and right 
of the line. The wavelength ranges we used for the maps are 
Bry: 2.16 - 2.17 yum; [Eelll] (2.242 yum): 2.2395 - 2.2438 yum; 
[Eelll] (2.3479 yum): 2.3460 - 2.3515 yum; [Eelll] (2.2178 yum): 
2.2154 - 2.2197 yum; CO-absorption: 2.2932 - 2.3015 yum. These 
wavelength ranges are several times larger than a typical radial 
velocity oflTset. This is true both for the stellar sources and the 
extended emission in the observed region, as the absolute values 
of radial velocity do not exceed 150 kms"^ (Paumard et al. 2004, 
2006; Zhao et al. 2009). The absorption line maps presented in 
Eig. [9] were created by extrapolating the continuum short- ward 
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of 2.29/im and subtracting it from the measured spectrum. The 
maps then were produced by plotting the negative value of the 
line integral over the strongest CO-bandhead between 2.290yum 
and 2.300yum. 

In Fig.[TOlwe show the radio image from Zhao et al. (2009) 
with line maps in contours over-plotted. The radio data and 
our imaging field spectroscopy data cubes have essentially been 
taken at the same epoch. In addition, we extracted spectra for 
specific regions as shown in Fig. O The radius of the extraction 
regions is 5 pixels (0.25"). In Fig. [TT] we show the integrated 
K^-band spectrum of IRS13N and for comparison the spectra of 
the nearby regions east and west of IRS13N that are dominated 
by diff'use emission. Integrated K^-band spectrum of the [Felll]- 
emission line blob north of IRS 13 N is presented in Fig. [121 As 
indicated by the shape of the [Felll] line at 2.3479yum this line 
might be a blend with the weaker Hell line that is close in wave- 
length. 

3. Results 

At a distance of 8 kpc the GC is well suited to derive proper 
motions of stars (Eckart & Genzel 1996, Ghez et al. 2000, 
Eckart et al. 2002) and dusty filaments (Muzic et al. 2007). In 
the framework of the GC, proper motions have mostly been 
used to derive a mass estimate of the potential the stars are 
moving in (e.g. Schodel, Merritt & Eckart 2009, Schodel et al. 
2005) or to explore the origin and dynamics of young and old 
stellar populations in the central cluster (Genzel et al. 2003, 
Levin & Beloborodov 2003, Ghez et al. 2005, Paumard et al. 
2006, Schodel et al. 2005, Muzic et al. 2008, Schodel, Merritt 
& Eckart 2009). However, the comparison of proper motions 
is also a unique tool to identify sources across diff'erent wave- 
length bands. In this paper we use this technique to identify stel- 
lar counterparts of dusty sources in the GC stellar cluster. 

3.1. Proper motions of the IRS13N and IRS13E sources 

To derive the K^-band proper motions of the IRS13N cluster 
sources it is necessary to work on images with a very high point 
source sensitivity for each epoch. In each epoch the faint K^- 
band counterparts of the significantly brighter L'-band sources 
had to be identified. A first identification for epoch 2002 was 
given by Eckart et al. (2004). The K^-band sources are faint com- 
pared to the IRS13E sources and deconvolution of point spread 
function (PSF) fitting is not suflftcient to provide source identifi- 
cations for all epochs. Therefore, we produced high-pass filtered 
images of the highest Strehl K^-band images. In these images 
the peak of the diff'raction limited core of the adaptive optics 
images is highlighted and the surrounding emission is strongly 
suppressed. 

The high-pass filtered image can be produced by calculating 
a smooth- subtracted image S as 

i: = I^G-I = (A-^E)^P^G-(A-^E)^P. (1) 

The symbol denotes the convolution operator, A denotes 
the distribution of stellar sources, and E represents the distribu- 
tion of extended emission. Here P is the point spread function 
and / the image that results from the convolution of P and the 
object distribution O = A-\- E. We use a narrow Gaussian G nor- 
malized to an integral power of unity. It is our experience that it 
is best to chose the Gaussian G to have a width that corresponds 
to about 30% to 50% of the full width at half maximum of a 



diff'raction limited core of a stellar adaptive optics image. The 
smooth- subtracted image can then be rewritten as 

i: = A>^(P>^G-P)-\-(E^P>^G-E^P). (2) 

Since the width of the extended emission is much larger than 
the width of the PSF (for the one-dimensional case see discus- 
sion on narrow dust features by Muzic et al. 2007), we can then 
infer that E^P^G - E^P and the high-pass filtered diff'erential 
image is 

S ~ A >^ (P >^ G - P) (3) 

(see also Sabha et al. 2010). If the image is clipped at zero 
flux density the remainder of P^G-P represents the top section 
of the difl'raction limited PSF core of the adaptive optics image 
and S is an approximation of an image of the point sources in the 
field. To improve the signal to noise for signals close to the spa- 
tial frequencies of the PSF we convolve the high passed filtered 
image again with the Gaussian G. 

This high-pass filtering method has the advantage that no 
estimate of the PSF is required. The PSF is typically variable 
across the imaging field and has a finite size. The correspond- 
ing apodization radius and residuals from stars that have to be 
removed from the PSF wings (via a median procedure or ex- 
plicit subtraction) are usually the cause for additional artifacts 
that spoil the image quality if high sensitivity is required. In 
addition, a deconvolution may suff'er from small number divi- 
sion problems at high spatial frequencies or it runs with a lim- 
ited number of iterations, or a loop gain and a stopping or con- 
vergence criterion. In case of PSF fitting a limited number of 
sources need to be identified. All of these restrictions and their 
potential drawbacks do not have to be considered using the high- 
pass filtering method. The noise property of the high-pass filter- 
ing is intimately linked to the original noise in the image. All 
noise contributions with spatial frequencies significantly lower 
than the characteristic cutoff' frequency of the filter (i.e. given 
by the FWHM of the Gaussian kernel) are strongly suppressed 
by subtracting the ffux conserved image from the original. Only 
(already existing) image power around and above the charac- 
teristic cutoff' frequency of the filter will remain. This includes 
power at the diffraction limit. Hence, there is no mechanism over 
which additional and spurious noise peaks can be introduced by 
the filtering process. High-pass filtering has, however, the dis- 
advantage that PSF imperfections at high spatial frequencies, as 
they may be introduced by the combined imaging optics (resid- 
ual astigmatism, trefoil etc.) are not corrected for. 
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Table 2. Relative and absolute coordinates of the lower left 
(south-eastern) corner of the images sections shown in the corre- 
sponding figure. The relative coordinates are given with respect 
to the SgrA* radio position in the first line. 



We derived the positions of individual sources in the high- 
pass filtered images from (1) Gauss-fits to the resulting source 
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profiles and (2) by determining the peak position on a sub-pixel 
level. Within the uncertainties (see below) both methods resulted 
in the same proper motion values. In order to provide an accurate 
relative positional reference system the low proper motion stars 
26-33 identified in Fig.[2land the May 2005 positions of stars El 
through E6 in the IRS13E cluster (Eritz et al. 2010) were used 
to establish a calibrated reference frame that allows us to com- 
pare the positions and proper motions to values in the literature: 
in particular Eritz et al. (2010) mainly for the IRS13E sources 
in K,-band, Muzic et al. (2008) for the IRS13E and IRS13N 
sources in L'-band, and Schodel et al. (2005) for the IRS13E 
sources and stars in the wider 3" diameter field around the clus- 
ter obtained in K^-band. 

The proper motions of sources a through k can be inspected 
in Eigs. |4l[5]and[6j Eor demonstration purposes the pixel scale 
of 27 mas has been decrease (by interpolation) by a factor of 
8. Grey scales have been adopted to the quality of the frame. 
Sources may in addition be influenced by blends with nearby 
background sources. In Tab. [3] we present the source name (col- 
umn 1), a reference (2), a source identification (3), coordinates 
relative to SgrA* and the corresponding uncertainties (4-7), and 
proper motion values and the uncertainties (8-11). The uncer- 
tainties for the positions and proper motion values were derived 
from the linear regression fit to the positions at all epochs. 

In Eig. [7] we compare the proper motion values of IRS13E 
and IRS13N sources obtained in L'-band (Muzic et al. 2008) to 
those obtained here in K^-band. We have compared our proper 
motion results with published values and quote the root mean 
square deviations in right ascension ctra. and declination (TDec. • 

These values have to be compared to the expected one- 
dimensional proper motion velocity dispersion value of about 
200 km/s at the location of the IRS13E and IRS13N clusters 
(Schodel, Merritt & Eckart 2009). 

Eor the IRS13E sources El to E6 we find good agreement in 
comparison to the results of Eritz et al. (2010) ctr a. = 11 km/s 
and (TDec. = 22 km/s as well as in comparison to the L'-band 
proper motion by Muzic et al. (2008) with ctra. = 38 km/s and 
CTDec. = 44 km/s. We regard the diff'erence with respect to Eritz 
et al. (2010) as a measure of our measurement uncertainty for 
bright sources in K^-band and the diff'erences to Muzic et al. 
(2008) a measure of the combination of our measurement uncer- 
tainties in K^-band and the corresponding L'-band uncertainties. 
The deviations in comparison to Schodel et al. (2005) are larger 
with ctr a. = 96 km/s and (TDec. = 77 km/s, but here the base- 
line in time used to derive the motions was only 3 years. Eor the 
IRS13N sources jStorjwc compare our K^-band proper motions 
to the L'-band values by Muzic et al. (2008) and find ctra. = 
65 km/s and (Tdcc. = 95 km/s. Correcting for the uncertainties 
derived for bright sources from the comparison to the data by 
Muzic et al. (2008) we obtain a residual uncertainty value be- 
tween the K^-band and L'-band proper motions of (Tirsi3n,r.a. = 
52 km/s and (TiRsuN,Dec. = 84 km/s that can be solely attributed 
to the faint co-moving IRS13N sources jStor] (see section |4l). 



3.2. Proper motions of the DSO 

In the central arcsecond the source density is highest in the over- 
all cluster. In this region we need highest angular resolution in 
addition to high point sources sensitivity to separate the faint 
and bright stars from each other and to search for a stellar coun- 
terpart of the dusty high velocity S-star cluster source reported 




0.66 (llight month) 

Fig. 4. K,-band identification of individual IRS13E and IRS13N 
sources at two different epochs and visualizations of the positions of 
the IRS13N sources a, j3, 6, 6, 77, and k at two different epochs. The 
source labels correspond to those in Tab. [3] 
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Fig. 5. Proper Motion of IRS13E and IRS13N sources. For source la- 
bels see |4] and [6] 



by Gillessen et al. (2012). In addition to the high pass filtering 
we therefore employed Lucy and linear Wiener filtering. A de- 
tailed comparison of the latter two algorithms in comparison to 
the CLEAN algorithm, which is sensitive to structures that are 
favored by high pass filtering, is given in Ott, Eckart & Genzel 
(1999). In the K^-band we used the same images as in the case 
of IRS13N. In the //-band we used the images employed for the 
HK-spectral index analysis of SgrA* presented by Bremer et al. 
(201 1) covering the years 2003 till 2008. In both bands we inves- 
tigated the images at the derived L'-band positions of the dusty 
S-cluster object relative to the corresponding positions of the star 
S2 and SgrA*. 

In addition we compared the data with detailed simulations 
of the field including important neighboring sources like S23, 
S54 and S63. For these sources the proper motions and in the 
case of S23 the curvatures, are listed in Tab.9 by Gillessen et 
al. (2009). From S63 we derived the proper motions from our 
H- and K^-band data, since the values listed by Gillessen et 
al. (2009) indicated a motion almost perpendicular to what we 
found. Our simulations clearly show that the DSO and star S63 
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0.^^'' (0.5 Ught months) 



Fig. 6. Visualization of the K^-band identifications of the IRS13N sources y, 6, 6, and ^ at 4 diff'erent epochs. The source labels correspond to 
those in Tab. [3] 



blend each other from 2002 until 2007. In 2007 the DSO starts 
to emerge from the S63 PSF and is clearly distinguishable from 
the neighboring stars. 

This results in K^-band identifications for all epochs since 
2007. We can use the 2008 H-band data to derive a lower limit 
for the DSO magnitude. Taking the background from the neigh- 
boring stars and the 3(T variation across the dififraction limited 
beam into account, we find a limit of m//>21.2. For a detailed 
identifica tion o f th e DS O in diflTerent epochs see Figs. lA.lllA.21 
IA.3[ IA.4[ IA.5I and IA.6[ In Fig. [141 we show the proper motion 
of the DSO projected on the sky in comparison to the orbital 
solution given by Gillessen et al. (2012). 



3.3. Proper motions in the wider field of the central star 
cluster 

The central parsec of the GC reveals diflTerent structures when 
observed at diflTerent wavelengths. In Fig. 14 shown by Eckart et 
al. (2006) these sources are denoted with D followed by a num- 
ber. In Fig. [T] (taken on July 6^^ 2004 in the L'-band) we show 
the proper motions of some of these sources in the central few 
arcseconds (see also Fig. lA.13l and. Fig. lA.14l in the Appendix). 
They are typically brighter than 14.6 in L'-band. Some of them 
appear to be circular symmetric on the sky (e.g. D2 and D5), 
others appear to be elongated (D3 and D4) or even filamentary 
as the source D8 at the eastern edge of Fig. [T] about -1.5" east 
and -0.5" south of SgrA*. The extended source D8 will be ex- 
cluded from our investigation since it is a very complex region 
in projection on the sky for which proper motions are difl&cult to 
obtain. Comparing Fig. 1 shown by Gillessen et al. (2009) to our 
Fig. [T] we can associate the position of the S-star S43 with the 
position of D2 (see also Tab. |4l). Within the uncertainties both 
objects also compare well in their Kg- and L'-band proper 
motions. For sources D3, D5, and D7 we could only determine 
L'-band proper motions. For D4 the data are given by Muzic et 
al. (2010). For D2, D6, S90, El, F2 the K,- and L'-band proper 
motions agree within the uncertainties given in Tab. IH 

The mosaics in L'-band were put into an astrometric frame 
(see Muzic et al. 2008). The H- and ^^-band data have been 
transformed into the coordinate system of the ^^-band mosaic 
obtained on 3P^ of July 2002, using the coordinate transforma- 
tion procedure described by Muzic et al. (2008). The positions 
of the 30 reference stars are evenly distributed over the field and 
were determined using the IDL code STARFINDER (Diolati et 
al. 2000) for the coordinate transformation. The known proper 
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Fig. 7. A comparison of K^- and L'-band proper motions in right as- 
cension and declination for the GC IRS13E and the IRS13N cluster 
(indicated by circles in dashed lines). The direction of the motion, i.e. 
S to N or E to W is labeled. The source labels correspond to those in 
Tab. [3] The K^-band motions derived here are labeled with the open 
dots. Eor IRS 13 and IRS13N the filled dots indicate the proper motions 
pubhshed by Eritz et al. (2010) and Muzic et al. (2008), respectively. 
The open circles depict the ^^-band proper motions described here in 
our paper. 



motions, to calculate the reference stars positions in the epoch 
of the reference frame, were taken from Schodel et al. (2009). 

The peak positions of the objects were determined using two 
dimensional Gaussian fits in each mosaic and band. The uncer- 
tainties in the graphs are the Icr uncertainty in the peak position 
and the uncertainties in velocities are the Icr-uncertainty of the 
linear fits. Figs. lA. 1 3l and lA. 14l show detailed maps of the imme- 
diate environment of the infrared excess sources close to SgrA* 
in the central stellar cluster. Details of the proper motions are 
shown in Figs. IA.15I and IA.16I and listed in Tabs. lA.li IA.2I and 
IA.3[ As an overview the results are summarized in Tab. (Hand in 
Fig. [TJ The figure shows that the motion in the diflTerent bands 
for each source can be related to each other, since the vectors 
point, within the error bars, in the same direction. 
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name 


ret 


ID 






Ao 


0-A6 










El 


smf 


1 


-2.955 


0.001 


-1.645 


0.001 


-143 


8 


-115 


5 


E2 


smf 


2 


-3.167 


0.001 


-1.730 


0.001 


-245 


8 


34 


5 


E3 


smf 


3 


-3.180 


0.001 


-1.532 


0.001 


-97 


18 


-14 


15 


E4 


sf 


4 


-3.206 


0.001 


-1.413 


0.001 


-235 


12 


-11 


4 


E5 


mf 


5 


-3.378 


0.002 


-1.502 


0.002 


-150 


32 


155 


28 


E6 


sf 


6 


-3.005 


0.001 


-1.420 


0.001 


-146 


9 


16 


14 


a 




8 


-2.707 


0.004 


-1.459 


0.003 


27 


63 


200 


41 






9 


-2.876 


0.002 


-1.255 


0.004 


-8 


31 


325 


58 


7 




22 


-3.140 


0.004 


-0.990 


0.003 


-208 


60 


262 


44 


5 




17 


-2.929 


0.004 


-0.825 


0.003 


-101 


63 


352 


41 


e 




13 


-2.889 


0.001 


-1.010 


0.001 


22 


14 


222 


20 


^ 




34 


-3.183 


0.004 


-0.812 


0.003 


-162 


65 


362 


53 




s 


19 


-3.134 


0.002 


-0.660 


0.003 


-114 


24 


194 


46 


K 


s 


21 


-3.094 


0.001 


-0.891 


0.002 


-172 


18 


-283 


28 






7 


-2.430 


0.001 


-1.697 


0.002 


31 


15 


126 


33 




s 


10 


-3.071 


0.001 


-1.235 


0.001 


-10 


18 


114 


22 






11 


-2.982 


0.004 


-1.162 


0.003 


-312 


56 


358 


54 




s 


12 


-2.836 


0.001 


-1.116 


0.001 


-180 


11 


-30 


22 




s 


14 


-2.757 


0.001 


-1.017 


0.003 


274 


17 


-261 


40 




s 


15 


-2.334 


0.001 


-0.898 


0.001 


-119 


19 


62 


19 




s 


16 


-2.486 


0.001 


-0.851 


0.002 


65 


19 


-249 


37 






18 


-2.955 


0.002 


-0.726 


0.003 


111 


27 


-269 


42 




s 


20 


-3.392 


0.001 


-0.746 


0.001 


-102 


9 


-130 


22 






23 


-3.253 


0.003 


-1.063 


0.001 


-155 


46 


105 


20 




s 


24 


-3.312 


0.001 


-1.188 


0.001 


22 


15 


-1 


20 




s 


25 


-3.544 


0.001 


-1.268 


0.001 


133 


15 


81 


20 






26 


-2.651 


0.004 


0.398 


0.002 


-436 


61 


-215 


30 






27 


-2.995 


0.002 


0.362 


0.002 


38 


24 


-201 


25 






Zo 


^ ^Q9 
-D.Dyz. 


U.UUJ 


U. JDZ 


n 009 


1 ZL9 


zLi 


1 1 s 

- 1 i J 


97 
z / 






29 


-4.21 \ 


0.002 


-1.314 


0.001 


81 


28 


-226 


21 






30 


-4.370 


0.001 


-1.638 


0.001 


-139 


20 


-71 


20 






31 


-4.049 


0.001 


-2.064 


0.001 


93 


22 


-145 


17 






32 


-1.607 


0.002 


-0.514 


0.001 


196 


24 


-57 


17 






33 


-1.346 


0.002 


-0.342 


0.001 


10 


25 


-270 


18 



Table 3. K^-band positions in arcseconds relative to SgrA* and proper motions in km/s of IRS13E and IRS13N sources at epoch 
2005.83. The abbreviations s,m,f stand for Schodel et al. (2005), Muzic et al. (2008) and Fritz et al. (2010), respectively. The 
identifications correspond to those in Fig.O 



name 


Aa 


+ 


AS 


+ 




+ 


V5 


+ 


H 


K. 


U 


H-K, 


K,-L' 


D2, S43 


-0.515 


0.010 


-0.114 


0.010 


+226 


44 


+480 


203 


18.7 


16.5 


13.5 


2.2 


3.0 


D6, S79 


-hO.699 


0.010 


-0.520 


0.010 


-112 


20 


+ 172 


50 


17.6 


15.4 


13.5 


2.2 


1.9 


S90 


-hO.585 


0.010 


-0.959 


0.010 


-127 


20 


+44 


60 


17.4 


15.3 


14.6 


2.1 


0.7 


Fl 


-0.321 


0.025 


-1.238 


0.025 


-265 


43 


+260 


50 


18.3 


16.0 


13.7 


2.3 


2.3 


F2 


-0.048 


0.025 


-1.546 


0.025 


-635 


200 


-77 


81 


18.0 


15.8 


12.9 


2.2 


2.9 


D3 


-0.623 


0.020 


-0.045 


0.020 


+410 


200 


+331 


134 


21.3 


19.8 


13.8 


1.5 


5.7 


D5 


-hO.547 


0.020 


-hO.030 


0.020 


+419 


150 


+468 


131 


19.8 


17.3 


13.4 


2.5 


3.9 


D4, S50, XI 


-0.540 


0.010 


-0.520 


0.010 


-52 


12 


+546 


15 


19.3 


17.1 


12.6 


2.2 


4.5 


D7 


-hl.338 


0.020 


-hi. 003 


0.020 


-132 


40 


-318 


40 


>23.4 


18.5 


14.1 


4.9 


4.4 


DSO 


-hO.187 


0.010 


-hO.065 


0.010 


-705 


170 


+560 


160 


>21.2 


18.9 


14.4 


>2.3 


4.5 



Table 4. Names, positional off'sets from SgrA* in arcseconds, proper motions in km/s of compact infrared excess sources in the 
central 2" field close to SgrA*. The data are given for epoch 2008. 



4. Discussion 

4.1. Identification and proper motions of IRS13N sources 

The derived stellar one-dimensional velocity dispersion at that 
location is 3 times larger than the root mean square diff'erence 
between the proper motions of the sources identified in the 
and L'-band. In addition, the diff'erence between our 2005.83 
one-dimensional IRS13N source positions in the K^-band and 
the L'-band positions by Muzic et al. (2008) at a mean epoch 
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Fig. 8. 2yL/m-band image (left; IB 2A5/um) and L'-band image (right; 
3.8yum) with the Bry emission (top) and CO absorption (bottom) line 
maps in contours over-plotted. 
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Fig. 9. 2yL/m-band image (left; IB 2.15 micron) and L'-band image 
(right; 3.8yum) with [Felll] -emission line maps in contours over-plotted 



closer to 2005 is of the order of 20 mas, i.e. two thirds of a K^- 
band pixel of 27 mas size. This is of the same order (or below) 
of what is expected as the maximum change in position due to 
proper motions of the sources during the time difference between 
these mean epochs. We conclude that over the past 10 epochs we 
used to analyze the motion of the IRS13N sources the K^-band 
sources can indeed be identified as being the co-moving coun- 
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Fig. 10. The 13mm VLA radio image from Zhao et al. (2009) with our 
Bry- and [Felll] -emission line maps in contours over-plotted. 



terparts of the L'-band sources. This also means that the L'-band 
sources are not unrelated (core-less) dust clouds as speculated 
by Fritz et al. (2010). 

It is likely that not all of the sources a to a: that had been 
discussed in Eckart et al. (2004) belong to the same association. 
Already in Eckart et al. (2004) we pointed out that k has sig- 
nificantly bluer colors than the other IRS13N objects. In fact it 
has not yet been identified with an L'-band source. Our K^-band 
proper motion analysis also shows that it is moving south and 
west rather than predominantly north as all the other IRS13N 
members. From this we conclude that k does not belong to the 
small IRS13N cluster of co-moving infrared excess sources. 

Source a is furthest south of the Eckart et al. (2004) sources 
and has also a velocity off'set from the other IRS13N sources. 
Here we give the L'-band velocities derived by Muzic et al. 
(2008) a higher weight, since in K^-band the source is appar- 
ently located in a crowded region as well as in the wings of the 
bright IRS13E member El. We conclude that a is also not likely 
one of the small IRS13N cluster of co-moving sources. We are 
therefore left with the sources j3 to r] that have been identified 
as co-moving and whose dynamics were analyzed in detail by 
Muzic et al. (2008). 

The source t] is also significantly brighter in //-band than all 
other IRS13N sources. If we assume that the IRS13N sources 
are embedded in their own dust disk or envelopes and that they 
all have a random orientation towards the observer this implies 
an average opening angle of 18076=30^ under which the source 
can be seen unobscured. From this one may conclude that the 
circum- stellar material almost entirely surrounds the objects. If 
the sources have disks it means that either not all of the material 
has settled into the disk, or that these disks are as thick as they 
are in diameter, or that they are significantly warped. All of this 
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Fig. 11. The integrated K^-band spectrum of 1RS13N (top panel), spec- 
trum of the nearby region Diffuse-1 west of IRS13N (middle panel) and 
Diffuse-2 east of IRS13N (bottom). The three spectra were normalized 
to the peak value of the Bry line in the 1RS13N spectrum. The spectra 
are continuum subtracted. 



Hel Felll Bry Hell Felll Felll 



Fig. 12. The integrated K^-band spectrum of the [Felll] -emission line 
source at the position of source K21 (Zhao et al. 2009) north of IRS 13 
N. This spectrum is continuum subtracted normalized to the peak of the 
Bry line. 



could be expected if the sources indeed belong to a dynamically 
young compact cluster of co-moving young stellar objects. 

The velocity dispersions of sources jStOT] corrected for their 
measurements uncertainties is of the order of 45 km/s for the L'- 
band data discussed by Muzic et al. (2008). For the less well de- 
termined motions derived from our K^-band identifications this 
value is about a factor of 1.5 higher. The implications for the 
possibly associated stellar mass of the IRS13N association has 
been discussed in detail by Muzic et al. (2008). The result is that 
if the IRS13N cluster is gravitationally bound then these veloc- 
ity dispersions imply a stellar mass that should be detectable. 




2008 




0,2 



arcsec 



Fig. 13. Discriminating the DSO against stars S57, S23, S54 and S63 
in 2008 in L'- and K^-band. Only an upper limit can be given in H-band. 



Hence, it is likely to be unstable implying that the association 
and its members are (dynamically) young. 
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0.3 0.2 0.1 0.3 0.2 0.1 

relative R.A. in arcseconds relative R.A. in arcseconds 



Fig. 14. Positions of the dusty S-cluster object relative to S2 and SgrA* 
plotted in comparison to the orbital tracks provided by Gillessen et al. 
(2012) (red for the DSO, blue for S2, SgrA* is at the position of the 
black dot at the origin). The results for the and L'-band are shown in 
panel a) and b), respectively (see data in Tab.[5l). Thin lines indicate the 
approximate expected positions of the DSO on its orbit from 2002.0 to 
2012.0 (based on our L'-band reductions and Gillessen et al. 2012). The 
epochs are color coded and the measurement uncertainty is shown as a 
black cross (see lA.Sl in the Appendix). 

4.2. Nature of the IRS13N sources 

In this paper we use the SINFONI line and continuum images 
for a qualitative interpretation of the nature of the extended and 
stellar sources only. A quantitative analysis of the line fluxes is 
beyond the scope of the paper and will be published elsewhere. 
The comparison between the Bry line emission and the near- 
infrared continuum emission of the stars in the IRS13E/IRS13N 
field in Fig. [8] shows that the line emission peaks on one of the 
brightest components (E2) of IRS13E, covers the entire IRS13E 
cluster and extends to the IRS13N complex. The strength of the 
emission is correlated with the presence of stars in the field and is 
clearly not only associated with the stars but comprises a diff'use 
component as well. The maps of the [Felll] -emission depicted 
in Fig. [9] show that the line flux is peaked both on IRS13E (E2) 
and IRS13N. The [Felll] flux show contributions from a difl'use 
component as well. In IRS13N the [Felll] integrated line flux 
density peaks to within <0.1" close to sources e and S. Zhao et 
al. (2009) find a radio source labeled as K23 that lies between 
6 and S and has to within 3-4cr the same velocity as the mean 
velocity of the NIR counterparts of e and S. 

The integrated [Felll] line emission is also peaked on a 
source about 0.15" north of t] with off'sets from SgrA* of 
-3.25"±0.03" and -0.45"±0.03" in right ascension and dec- 
lination, respectively. This position agrees with that of source 
K21 identified by Zhao et al. (2009) to within about 0.1". The 
authors find a proper motion velocity of Vq,=-224±16 km s"^ 
and v^=-421±20 km s"^ towards the NNW. In the L'-band we 
do not find a well defined point source at the position of K21. 
The continuum flux density extension from IRS13N at that posi- 
tion is consistent with a source that is about a magnitude fainter 
(i.e. L'~12) than source t]. At K^-band no clear association to a 
source brighter than K~16.3 can be established. This indicates 
that K21 is either associated with a stellar source even more red- 
dened than the IRS13N sources or that the source K21 is a core- 
less dust source bright in [Felll] -emission. There is also exten- 
sion in the integrated line emission of Bry covering the position 
of K21 . A spectrum of the infrared flux density at the position of 
K21 is shown in Fig.[T2l 

The comparison between the CO absorption and the near- 
infrared continuum emission in Fig. [8] shows that the IRS13N 
complex is clearly located at a minimum of the CO-absorption 
in the IRS13E/IRS13N field. The absorption is concentrated on 



a few stars in the field to objects about 0.5" east, west and north 
of IRS13N. The minimum in the absorption strength centered on 
IRS13N comes close to the background value due to the overall 
central stellar cluster. Given the contribution through PSF wings 
of these stars at the location of IRS13N this is fully consistent 
with the statement that most of the IRS13N sources (about 0.5" 
north of IRS13E) do not show any significant CO line absorp- 
tion. 

The comparison of the radio data by Zhao et al. (2009) 
(Fig. [T0|) shows that the radio continuum emission peaks to 
within <0.10" at the peaks or contour line excursions of the 
[Felll] line emission. The radio continuum is more peaked on the 
IRS13E source El. The correspondence to the [Felll] -emission 
(in particular the 2.242yum line emission) is best. The Bry emis- 
sion is clearly more extended than the bright radio continuum 
emission, highlighting that the radio emission is more dominated 
by the compact stellar sources than the more diff'use emission 
traced by Bry. 

The combined spectra on IRS13N and the two off'-positions 
east and west (Figs. [3] and [TT]) demonstrate that the emission in 
the Bry, H2 and in particular in the [Felll] lines is brighter on 
IRS13E than in the off'-positions. The spectra also demonstrate 
that the integrated CO-band head absorption towards IRS13E 
is as strong as towards the off'-positions. This indicates that the 
CO-absorption at this position can be fully explained by the CO- 
absorption in the PSF wings of the bright nearby late type stars 
and the overall stellar cluster background. 

4.3. NIR colors and spectra of the IRS13N cluster 

Compared to the immediate surroundings (except for the 
IRS13E cluster) the sources in the IRS13N region are clearly 
correlated with a minimum in CO absorption and an excess in 
Fe- and H-recombination line emission. While the excess in 
emission may in part be due to a local concentration of gas 
associated with the 13N association or possibly with the mini- 
spiral, the lack of CO absorption can be attributed to the pre- 
dominantly early spectral type as suggested from the NIR colors 
of the sources (Eckart et al. 2004). 

While near-infrared CO absorption and emission is often 
present, it is not a necessary feature to be observed towards 
young massive stars. Out of 201 young massive stars in M17 
Hoff'meister et al. (2006) find 30% to have featureless spectra. 
Half of those show an infrared excess which corresponds to ba- 
sically the same percentage of the remaining 60% of the sources 
that show CO in absorption or emission. Almost a third of the 
featureless sources also tend to show X-ray emission and a fifth 
shows simultaneous infrared excess and X-ray emission. 

Featureless sources have the tendency to be brighter in the 
near-infrared compared to other young stellar objects. From the 
distribution in Fig. 3 in Hoff'meister et al. (2006) we can con- 
clude that featureless objects are typically up to 2 magnitudes 
more luminous compared to CO emission or absorption sources. 
They are apparently later than spectral class B3. 

Compared to previous CO bandhead studies in star-forming 
regions, the investigation by Hoff'meister et al. (2006) is the most 
extensive and a rather unbiased one. Casali & Eiroa (1996) argue 
that Class II sources tend to show CO absorption, while Class I 
sources are featureless. The thermal continuum emission from 
dusty in-falling envelopes can exceed the combined stellar and 
disk photospheric emission by almost an order of magnitude 
(Calvet et al. 1997, Hoff'meister et al. 2006). Hence, this pro- 
cess can result in featureless spectra. In addition, the veiling of 
an envelope will also weaken the CO absorption lines. The pres- 
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ence of envelopes with large scale heights is consistent with the 
presence of a single bright source 7/ in addition to 5 faint objects 
as elaborated above. The presence of a disk, on the other hand, 
is likely to amplify the CO feature. Therefore, Hoffmeister et al. 
(2006) conclude that the 60 featureless objects in their sample of 
201 sources in Ml 7 are very likely Class I sources which have 
just started to build up an infrared excess that is progressively 
moving towards the MIR and NIR domain. 

From model calculations Hoffmeister et al. (2006) infer that 
the observed CO absorption is most likely a sign of heavily 
accreting proto-stars. The mass accretion rates may be above 
10"^ Moyr"^. High accretion can be provided through angular 
momentum loss within a disk. We conclude that the absence 
of that feature may indicate that these stars are currently not 
strongly accreting or have not yet reached that phase. This would 
be the case if the luminosity of the in-falling envelope is strong 
against the emission from the disk. 

We infer from the combination of the fact that the stellar 
sources are in a co-moving group (Muzic et al. 2008, Eckart et 
al. 2004 and the proper motions and the spectroscopy presented 
here - see below) that the IRS13N sources may be featureless 
young stellar objects, like Class I sources, in which a bright NIR 
disk has not yet been formed but in which the NIR/MIR excess 
is still dominated from surrounding and in-falling material. 

4.4. The age of the IRS13N sources 

Possible stellar identifications of the IRS13N sources have al- 
ready been discussed extensively by Eckart et al. (2004) and 
Muzic et al (2008). The fact that our data supports the finding 
that the sources (J3 to 77; see below) are indeed co-moving poses 
serious problems to source identifications that involve stars with 
ages >3 Myrs. Here we highlight a few more aspects concerning 
the identification of the sources with young stars: 

Broad band colors can give information on the nature of 
sources. In Fig. [16] (left) we compare the colors of dusty objects 
in the central few arcseconds of the Milky Way to the colors of 
young stellar objects as given by Hoff'meister et al. (2006). The 
sources of IRS 1 3N are higher extincted than the remaining dusty 
objects in the field. Their dereddened colors are in agreement 
with those young dust enshrouded objects (details in Eckart et al. 
2004). The assumption (see section |43l) that the IRS13N cluster 
members are Class I sources implies that they are young. In a 
recent study of stellar and circum- stellar properties of six Class 
I proto-stars Prato et al. (2009) estimate their ages as <2 Myr. 
However, even the ages of young Class Il/Class III sources are 
estimated to range between 0. 1 and 1 Myr (Green & Meyer 1995, 
Luhman & Rieke 1999, see also van Kempen et al. 2009). This 
is consistent with the assumption of Eckart et al. (2004) that the 
IRS13N sources are good candidates for being YSOs or young 
Herbig Ae/Be stars with ages that may range between 0.1 and 
about ~3 Myrs. Molinari et al. (2008) find that the time scales 
of the formation of stars is of the order of 1-4x10^ years, with 
shorter times for higher masses (6 to 40 M©) (Hillenbrand et al. 
1992, Fuente et al. 2002, Ishii et al. 1998). The orbital time scale 
for the combined motion of the entire IRS 13E/IRS13N com- 
plex and the mini-spiral gas of the order of -10"^ years (see also 
Muzic et al 2008) is small compared to the plausible ages of the 
massive young stars. We also assume that their associated en- 
velopes and disks can survive the hostile GC environment for at 
least 10"^ years (see discussion in Scally & Clarke 2001). 
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Fig. 15. Decomposition of the DSO spectrum including our K^-band 
detection and H-band limit. Here we demonstrate that a mixture of dust 
and stellar contribution is possible. The points correspond to the L- and 
K5 -band magnitudes, and H-band upper limit from this work, and the 
M-band measurement and K-band upper limit of Gillessen et al. (2012). 
Red and blue dashed curves show also their 550 K and 650 K warm 
dust fits. In solid blue, red and magenta lines the emission from three 
different possible stellar types of the DSO core are plotted. Any of these 
stars embedded in 450 K dust (solid brown line) can produce the black 
line which fits all the NIR DSO photometric measurements. 



4.5. Broad band spectrum of IRS13N 

Inspection of the 1.2" angular resolution 1.3mm map of the 
mini-spiral presented by Kunneriath et al. (2011) show that the 
IRS13N complex is associated with 1.3mm continuum emis- 
sion. Just north of IRS13E it shows itself as a compact source 
component with a flux density of 20+3 mJy (see Fig. IA.9I ). 
Arithmetically de-convolved with the beam its size is less than 
0.6 arcsec. In NE direction it is partially blend with IRS13E and 
a northern neighboring mm- source. In the 13mm map presented 
by Zhao & Goss (1998; Fig.3/PLATES therein) we can read off 
the map in their Fig. 3 an extended flux density of 1 1 ±2mJy with 
a compact source of about 3mJy. These flux densities are con- 
sistent with a dominant non-thermal stellar wind contribution 
with an approximate mass load of 10"^ M© yr"^ (e.g. Monies 
et al. 2009; Lang et al. 2005). Such a mass load is typical for 
young luminous stars (Drake & Linsky 1989, Kudritzki et al. 
1999, Crowther 2001). 

Eckart et al. (2004) quote H, and L' flux densities for 
individual IRS13N sources. Viehmann et al. (2005, 2006) give 
NIR/MIR flux density values obtained in an ^^1" diameter cir- 
cular aperture on IRS13N. In the J- and //-band and probably 
even in the K^-band these 1" aperture flux densities are domi- 
nated by confusing fore- and background sources of the central 
cluster and the values by Viehmann et al. (2006) can probably 
only be taken as an upper limit. Due to the very red colors of 
the sources the 1" flux densities at longer wavelengths are most 
likely characteristic for the IRS13N sources. However, at these 
wavelengths dust emission becomes dominant. Our investigation 
of the K^- and L'-band proper motions shows that the dust emis- 
sion dominated L'-band flux densities can be fully associated 
with the stellar sources found at shorter wavelength. Hence, we 
assume that a dominant portion of the emission longward of the 
/T^-band can in fact be attributed to the IRS13N sources. 
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Table 5. The K^, and L'-band positions of the dusty S-cluster object (DSO) and the star S2 relative to the SgrA* position. Typical 
uncertainties range between 13 mas and 25 mas for the DSO and (see Fig. IA.8l) less than 10 mas for the position of S2. 
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2004.00 
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0.000 


0.000 


DSO 


0.30 


2002.00 
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1.36 


2006.59 
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Table 6. Epochs, flux ratios, pos itio ns, pr oper motions and accelerations used to simulate the DSO and the neighboring stars with 
results shown in Eigs. lA.ll IA.2l and IA.3I The flux ratios have been calculated with respect to the mean flux derived (m^=17.6) 
from the K^-band magnitudes of S23, S54, S57, and S63 as given by Gillessen et al. (1999). They correspond (with an estimated 
uncertainty of 30%) to the relative fluxes suggested by the modeling presented in the figures listed above. Erom the quality of the fit 
in these figures we estimate that the uncertainties are 10 mas for the coordinates, 1 mas yr"^ for the proper motions and 0.1 mas yr"^ 
for the accelerations. 



4.6. Proper motions and nature of the central dusty sources 

Examination of the proper motions of the dusty sources in the 
central few arcseconds shows that 5 of the 10 L'-band excess 
sources identified in Eckart et al. (2004, see Eig.[T]) have K^-band 
counterparts. The color-color diagram in Fig. M (left) shows 
that the colors of dusty objects in the central few arcseconds 
of the Milky Way are bluer than those of the IRS13N sources. 
Correcting for the Galactic foreground extinction of approxi- 
mately 27 magnitudes in the visible places them at the location 
of the young stars discussed by Hoflfmeister et al. (2006). The 
IRS13N sources are higher extincted than the remaining dusty 
objects in the field. Correcting for the Galactic foreground ex- 
tinction places them at the location of extincted young stellar 
objects in Eig [16} (left). Sources /3 to r/ in IRS13N are there- 
fore excellent candidates for being young stars with an average 
extinction around Ay ~ 27. 

4. 7. Pure dust versus photospheric emission 

Here we investigate if the NIR emission of the infrared excess 
sources can be due to pure dust sources or if significant contribu- 
tions from photospherice emission of stars need to be included. 
If the objects were pure dust sources then their observed colors 
must be consistent with reddened colors of pure dust sources. 
Eor the IRS13N this apparently not the case. The DSO is clearly 
detected in the K^-band at m^=18.9. In H-band we find an up- 
per limit m//>21.2 if we consider the background level from the 



neighbouring stars S23, S54 and S63 in 2008, and the 3cr vari- 
ation at the position of the DSO in K^- and L'-band. With the 
DSO L'-band brightness of m^' = 14.4+0.3 obtained from the ob- 
serving epochs listed in Tab. [T] we get H-K^=2.3 as a lower limit 
(Eig. [T6l right). This results in a measured color K^-L'=4.5. With 
A//=4.21, A/^^=2.22 and Al' = 1.09, Am~Al' (e.g. Viehmann et 
al. 2005, Lutz et al. 1996) the absolute magnitudes can be found 
and plotted in a spectrum (Eig.[T5l). In Eig. [T6l (right) we com- 
pare the data to the black body line and extincted colors of dif- 
ferent mixtures between a stellar and a dust contribution (Glass 
& Moorwood 1985). The pure 550 K dust colors (K,-L'~3.0) 
reddened by Ay =27 are consistent with the K^-L' color of the 
DSO. However, (with the current H-K^ color limit and with po- 
tential H-K^ colors that range up to the pure 550 K color) the 
dereddened positions of the DSO show that the DSO colors are 
also consistent with a mixture between a stellar and a dust con- 
tribution. Withn its current (H-K^ color limited) position in then 
diagram, mixtures with 60% stellar emission and 40% dust with 
temperatures less than 550 K are possible. If the DSO will be 
detected half way between its current H-K^ and that of pure 
550 K dust (striped circle in Eig.[T6]), 20% stellar and 80% 550 K 
dust contribution are possible. In Eig. [15] we show a possible 
spectral decomposition of the DSO using the M-band measure- 
ment by Gillessen et al. (2012). In addition to their 550 K and 
650 K fit and their m^^ = 17.8 limit (dashed lines) we plot our K^- 
band DSO detection and the H-band limit. The integrated flux 
amounts up to 30 L© . We show a decomposition in which we 
assume that 50% of the current band flux is due to a late 
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Fig. 16. Left: Colors of dusty objects in the central parsec of the Milky Way to the colors of young stellar objects as given by Hoffmeister et al. 
(2006). The label COAS and COES indicate that the sources show CO absorption or emission in their spectra, respectively. Right: Colors of dusty 
objects in the central parsec of the Milky Way compared to colors of pure dust (temperatures given in red numbers) and mixtures of a stellar and 
a dust contribution (fractional numbers given in black; see also Glass & Moorwood 1985). Here we demonstrate that a mixture of dust and stellar 
contribution is possible. Because the H-K5 color of DSO is an upper limit, its NIR emission can be explained as pure dust at 550 K reddened by 
~ 27 mag, or as a star embedded in warm dust, e.g., if its color were detected to be H - Ks ~ 3.5 mag (dashed circle), after dereddening by 27 mag, 
its emission would be consistent with 20%stellar and 80% dust contributions. 



dwarf (blue line) or an AF giant or AGB star (magenta line). We 
use dust at a temperature of 450 K plotted in red (see deredden- 
ing of the DSO in the color color diagram in Fig. [16] right). The 
sum of the spectra is shown by the thick black line (dashed at 
short wavelengths for the AF giant/ AGB case). 

Similarly, after dereddening the IRS13N sources for 
Ay=21^ they miss the pure black body line by about 2 mag- 
nitudes (see Fig. [16] left). For fully dereddened objects an offset 
of 2 magnitudes in K^-L' calls for temperatures of about 800 K 
or below with flux density contributions of a pure dust compo- 
nent of less than 50% (see also Fig. 7 by Glass & Moorwood 
1985). This temperature estimate is consistent with the overall 
IRS13N temperature around 650 K given by Fritz et al. (2010). 
This fact clearly implies that the IRS13N sources cannot be pure 
dust emitters but must have a stellar core. This is also fully con- 
sistent with the fact that the colors of these objects are similar to 
those of early dust enshrouded stars with colors shown the dia- 
gram in Fig. [16] (left). Hence, for IRS13N sources the evidence 
for photospheric emission from a star is strongly supported by 



the H-band detections (Eckart et al. 2004 Tab. 2 - here all sources 
except 6 and y could be identified). 

Inspection of Fig. [16] (left), the black-body line shown in 
Fig. [16] (right) and comparison to Fig. 7 by Glass & Moorwood 
(1985) shows that within the uncertainties the sources Fl and 
D7 are good candidates for being pure dust sources. Sources D6 
and S90 will have mostly stellar colors after dereddening with 
a possible contribution of dust >1000 K. Sources F2, D2, D3, 
D4, D5 are clear candidates for being objects that show a mix- 
ture from photospheric emission and emission from dust with 
temperatures well below 1000 K. 

4.8. Statistical robustness of tlie DSO identifications 

Here we analyze the statistical robustness of the DSO identifica- 
tions in the H- and -bands. In Sabha et al. (2012) and Eckart 
et al. (2012) we tabulate the probability for finding a blend star 
above the confusion limit in a single K^-band spatial resolution 
element in the overall region of the S-star cluster and at the po- 
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Fig. 17. Comparison of DSO images in the a) Ks- (13.2 mas pixel scale) and b) L'-band (27 mas pixel scale). The image section size is T'xl". 
In c) and d) a Gaussian of 98 mas FWHM and a PSF from the image have been scaled to the DSO peak brightness and subtracted at its position, 
respectively. An astrometric grid of red circles marks some bright stars that can be used as positional references to compare the images at different 
wavelengths. 



sition of SgrA*. With two independent methods we determined 
the central K^-band confusion limit observationally. Sabha et al. 
(2010, 2012) employ successive subtraction of stars in the S-star 
cluster. Witzel et al. (2012) use the shape of the SgrA* flare am- 
plitude histogram. The otherwise straight power spectrum with 
which this histogram can be described turns over at a flare flux 
density at which a mere flare detection becomes a flare flux den- 
sity measurement (details in Witzel et al. 2012). This indicates a 
confusion limit in the range of 0.5 to 1.5 mJy. 

The result of the statistical investigation of the S-star clus- 
ter by Sabha et al. (2012) is that for a range of observationally 
supported parameters that describe the structure and luminos- 
ity function of the cluster, the authors can give the probabil- 
ity of finding a blend star in a single spatial resolution element 
above the confusion limit. Here we use a power law index of the 
projected spatial distribution of Qf=0. 30+0.05 and a K-band lu- 



minosity function exponent of KLF=0. 18+0.07. Projecting the 
S-star cluster onto the sky blend stars are created by line of 
sight overlaps of fainter stars within a finite angular resolution 
element. The probability of finding blend stars in the central 
1.3x1.3 arcsec^ overall can be as high as 1.0 and reaches a value 
of about 0.5 at the position of SgrA* itself. In the following we 
use these values although they can be much lower if the imme- 
diate vicinity of SgrA* is emptied from low mass stars through 
stellar dynamical processes. 

The probability of finding a K^-band source in an L'-band 
resolution element (that corresponds to 2.8 K^-band resolution 
elements) at a specific predicted position within the S-cluster is 
about 170 times smaller that unity, i.e. 6x10"^. As explained in 
Sabha et al. (2012), the livetime of these blend stars (also in- 
dicated through observations) is of the order of 3 years during 
which they dissolve again into individual objects with flux densi- 
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ties below the confusion limit due to their proper motions. If one 
wants to explain the 2007-2012 K^-band identifications of the 
DSO by blend stars one requires about 2-3 of these transient ap- 
parent objects. This happens with a probability of 3x10"^. Since 
there is no guarantee that apparent proper motions during the 
dissolving phase of the blend stars follow the observed orbit and 
since we have chosen upper limits to start with these low proba- 
bility values can be regarded as a safe upper limits for a scenario 
in which the DSO can be explained as the result of blend stars. 
Even if a blend star is assumed, one can think of at least 5-10 
proper motion velocity values and 5-10 directions that would be 
expected given the accuracy of 13-25 mas and the involvement 
of 2 to 6 epochs. This will lower the probabilities quoted above 
by a factor of 25-100. This implies that in a combined probabil- 
ity of being confused by a comoving blend star less than 10""^. 

However, the number of bright, single stars in the central 
region is high. We can therefore consider the likelihood of find- 
ing faint stars above the confusion limit at the position of the 
DSO over a time span of several years. Erom the KLE analysis 
by Sabha et al. (2010, 2012) we infer a maximum of about 40 
stars fainter than 16 magnitudes (dereddened) within the central 
1.3"xl.3". This implies 0.2 sources per resolution element and 
observing epoch. Eor 6 consecutive epochs this results in a prob- 
ability of the order of a percent to find a random configuration 
of the central cluster such that a faint star is found at consecu- 
tive positions given by the L'-band detections. If we include the 
immediate neighborhood (i.e. 9 to 25 surrounding resolution el- 
ements) the probability can rise up to 25%. With a limited accu- 
racy of determining the direction and value of the proper motion, 
this implies that the likelihood of finding a star that is close in 
position and velocity is in the percentage range. This may ex- 
plain the confusion of the DSO with S63 in the interval 2002 to 
2007. 

In summary, the contamination with a bright field star over 
the entire time interval 2002 to 2012 is likely (around a few per- 
cent) but being confused by blend stars is very unlikely (<10"^). 
Eor the IRS13N and all other inferred excess sources discussed 
here the projected stellar densities are even lower than those 
close to the center, such that serendipitous identifications are 
also not very likely. 

4.9. Proper motions and identification of ttie DSO 

The analysis presented here shows that the dusty S-cluster ob- 
ject reported by Gillessen et al. (2012) does indeed have a K^- 
band counterpart. Its H-K^ -colors indicate that it cannot be ex- 
cluded that the DSO is an embedded star rather than a pure dust 
cloud. Gillessen et al. (2012) showed that both in proper mo- 
tions and radial velocities the trajectory of the source perfectly 
fits a closed elliptical orbit. The authors interpret this source as 
a dust cloud which is approaching SgrA* and will be disrupted 
during peri-bothron passage. The combination of the fragility of 
the cloud and the fact that it is on a bound elliptical orbit implies 
that the object must have been on that orbit for a couple of rev- 
olutions. In the harsh environment of the central stellar cluster 
the evaporation time scale for a pure dust cloud is incompati- 
ble with this fact. The dust cloud would disappear before it can 
establish a relaxed orbit. If the object is on a bound orbit as a 
cloud and since it will be heavily stretched out by tidal forces 
during peri-bothron passages it should not be as compact as it 
currently appears. Variations in the shape should (with lower 
angular resolution than in the H- and K^-band) be traceable in 
the bright L'-band continuum where it appears to be rather com- 
pact (Gillessen et al. 2012). In the K^-band the DSO is close to 



the confusion limit (see discussion in subsection 14. 81) . Given the 
strong and extended Bry recombination line emission from the 
mini-spiral the situation may be similar for the DSO Bry line 
emission. As in the case of the continuum identifications this 
may influence the brightness, position, shape and velocity field 
of the source as seen in the light of the line emission (Gillessen 
et al. 2012). 

If the DSO is indeed a dusty star, then it may develop a 
bow-shock if it passes through an accretion wind from SgrA* 
quite similar to the sources X3 and X7 described by Muzic et al. 
(2010). Since the DSO as a dusty source is an obvious probe for 
strong winds possibly associated with SgrA*. Its presence and 
the fact that it has not yet developed a bow- shock, although it 
is already closer than X3 and X7, indicates that the wind from 
SgrA* is highly un-isotropic and possibly directed towards the 
mini-cavity (see Scenario IV below and Muzic et al. 2010 for 
further discussion). 

Comparing the luminosity and the magnitude of the DSO 
(Tab. |4]) we obtained in the H- and K^-band to the absolute 
magnitudes of OB and Be stars as presented by Wegner (2006) 
and stars of luminosity class I-V, we find that the DSO could 
preferentially be a K-dwarf (V) or a E-dwarf/sub-giant (V-III) 
at no additional extinction intrinsic to the surrounding dust 
cloud or an A-dwarf/sub-giant (V-III) at about 30 magnitudes 
of intrinsic visual extinction (see Eig. O. The AGB or early 
evolutionary phase may be included here in order to support the 
dusty envelope. These identifications point towards a slightly 
later type as compared to the earlier type BO-2.5 V stars in 
the S-cluster (Martins et al. 2008). This would imply a stellar 
mass of a few solar masses. However, the combination of local 
extinction (including geometrical eff'ects like shadowing and 
volume filling factors of the extincting material) and scattering 
may alter that identification and it could be even earlier or 
more luminous as well. This leads us to discuss four additional 
scenarios: 

Scenario I: Since the DSO is prominent in its hydrogen 
line emission (Gillessen et al. 2012) with a Bry line with of 
100 km/s in 2008 but rather faint in the H and continuum 
bands (see Tab. |4]) it may very well be a dusty, faint narrow 
line Wolf Rayet star. They are known to be less luminous than 
most other Wolf Rayet or OB giants and there are a few of these 
objects in the GC field (see e.g. Moultaka et al. 2005, Paumard et 
al. 2006). The Wolf Rayet stars IRS7W with m^=13.1, IRS7SE 
with m^=13.0 (Martins et al. 2007) and WR2 with m^=12.9 
(Moultaka et al. 2005) are amongst the faintest so far known in 
the Galactic Center cluster (for faint WR stars in the Milky Way 
see also Shara et al. 1999, 2012, Mauerhan, van Dyk, Morris 
201 1 with survey limits around m/^=15.0). Given those apparent 
magnitudes the DSO would have to be an exceptionally faint 
WR star. The earliest WR stars with dense winds and strong 
free-free excesses may belong to the intrinsically faintest of the 
WRs (Mauerhan, Van Dyk, Morris 2011). As indicated by the 
increasing Bry (Gillessen et al. 2012) the harsh environment 
of the S-star cluster and the possibly repeated interaction 
with the SgrA* black hole may have altered the stars mass 
loss or shell size and hence its luminosity in the thermal infrared. 

Scenario II: Dong, Wang & Morris (2012) have carried 
out a multi-wavelength study of evolved massive stars in the 
GC. In their Eig. 12 they show a Paa equivalent line width (EW) 
vs. K^-[3.6yum] plot of these objects. Erom this plot it becomes 
obvious that for WC stars and the fainter WNL stars/OB super- 
giants the EW values are smaller than lOOA for K^-[3.6yum]>2 
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magnitudes and even drop significantly with further increasing 
infrared excess reaching values of up to K^-[3.6yum]~5. Dong, 
Wang & Morris (2012) suggest that the free-free emission from 
the strong winds of the WNL stars/OB giants could potentially 
dominate the emission in the mid-infrared (Wright & Barlow 
1975). This could explain the positive EW vs. K^-[3.6yum] 
correlation they find. These properties fit well with the DSO if 
we assume that K^-[3.6yum]^^K^-L' and that a small EW in the 
Paof line is correlated with a small EW of the Bry line. The Bry 
line width in 2008 indicated in Fig. 2 by Gillessen et al. (2012) 
is of the order or smaller than 100 km/s in 2008 and 200 km/s 
in 2011. Taking the Galactic foreground extinction into account 
its position in the color-color diagram in Fig. [161 (left) comes to 
lie in the WNL stars/OB super-giants domain with small EW 
and high infrared excess as described by Dong, Wang & Morris 
(2012). 

Scenario III: Murray-Clay & Loeb (2012) propose that 
the dusty S -cluster object is a proto-planetary disk that has been 
brought in from a young stellar ring. However, to change the 
ellipticity of the previously rather circular orbit in one single 
event to the high ellipticity orbit on which it is now, requires 
a similarly violent interaction as the upcoming peri-bothron 
passage. If the object would be only a core-less dust cloud 
this event would therefore lead to an early destruction of the 
object at the time it enters the current orbit. If the source is 
indeed a core-less dust cloud the only alternative is that the 
ellipticity of its orbit has been changed gradually. However, this 
implies many peri-bothron passages, imposes a conflict with 
evaporation timescales and will most likely lead to destruction 
as well. The only alternative is indeed that the dusty S -cluster 
object has been formed on an orbit rather similar to the current 
one or that it has been brought in gradually from outside the 
S-star cluster as proposed earlier by Murray-Clay & Loeb 
(2012). This hypothesis is also supported by the analysis in the 
color-color diagram shown in Fig. [161 (left). They show that 
there are two objects, D3 and the DSO, which are dominated by 
their dust emission in the L'-band as it would be the case if they 
were relatively young stellar objects. 

Scenario IV: The DSO may be similar to the bow-shock 
sources X3 and X7 (Muzic et al. 2010) but has not yet devel- 
oped an obvious bow- shock structure. For X3 and X7 Muzic et 
al. (2010) discusses possible stellar counterparts including the 
possibility of late type and main sequence sources. However, 
the preference is clearly on late B-type main sequence stars (B7- 
8V) and low-luminosity Wolf-Rayet (WR) stars as well. Other 
explanations have problems: The central star of a planetary 
nebulae will remain at the required brightness for less than 1000 
yrs before it becomes too faint, and a main sequence star cannot 
be an appropriate source of a dust-rich envelope as apparently 
observed for the DSO in the L -band. 

4.10. Will the DSO be disrupted? 

Gillessen et al. (2011) report that the Bry emission of the DSO 
is spread over 200 mas. Similar sizes are also predicted in the 
models presented by Burkert et al. (2012). In Fig[T71we com- 
pare images of the DSO in the Ks- and L'-band. We also show 
the L'-band images after removal of the DSO using a Gaussian 
of 98 mas FWHM (i.e. the diflraction limit in the L'-band) and 
a star located 0.5" west and 0.42" north of SgrA* (i.e. a PSF 
obtained from an image section close to SgrA*; see Fig.[T]). The 
analysis of the images shows that >90% of the DSO emission 



at 3.8yum wavelength is compact (FWHM<20 mas). Only up to 
10% of the flux density contained in the compact part of the DSO 
can be extended on the scale size of the PSF. This indicates that 
the warm (-550 K) dust emission is very compact in compar- 
ison to the hotter (-10"^ K) hydrogen recombination line emis- 
sion. Alternatively the MIR continuum emission contains a sig- 
nificant compact free-free contribution (see below). In this case, 
however, one would expect that the recombination line emission 
is dominated by contributions from this free-free component and 
should hence be more compact than observed. 

Since the K^-band identifications of the DSO suggest that it 
can also be associated with a star it is only its dusty envelope 
that can potentially be disrupted. In order to investigate that pos- 
sibility we have to determine the location of the Lagrange point 
LI through which mass can be transferred and compare it to the 
size of the dust shell or disk. In Fig. [181 we show a sketch demon- 
strating the location of the DSO and the motion of the Lagrange 
point LI with respect to it. In panel a) in Fig. [181 we show the 
situation about 1 year before the peri-bothron passage. The LI 
point is rapidly approaching the DSO due to its orbital motion to- 
wards SgrA*. At peri-bothron LI will also be closest to the star 
(panel b) in Fig. [181). If the star has a mass of about IM© the 
separation of LI from it will be about 0.1 AU. For a Herbig 
Ae/Be stars with 2-8Mo that distance will be 0.2 and 0.5 AU. 
For a typical S-cluster stellar mass of ~20-30Mo the separation 
will be closer to an AU. Interferometrically determined typical 
inner ring sizes for young Herbig Ae/Be and T Tauri stars can 
be as smaU as 0.1-1 AU (Monnier & MiUan-Gabet 2002), how- 
ever, a disk or shell may be much larger than this. Based on its 
MIR-luminosity Gillessen et al. (2012) determine a source size 
of about 1 AU. Hence, during peri-bothron passage a significant 
amount of the dusty circum- stellar material may get beyond LI 
and start to move into the Roche lobe associated with SgrA*. 

However, this peri-bothron passage will only last for a few 
months after which the DSO will move away from SgrA* at a 
speed around a 1000 km/s within the first year. After one year 
the separation of the LI -point from the DSO is already about 
1000 times larger than during the peri-bothron passage (panel c) 
in Fig.[T8l). This implies that any cold dusty material that passed 
the LI towards the SgrA* will be overtaken again by the reced- 
ing LI -point. However, the time will not be sufficient to com- 
pletely disrupt the dust shell or disk of the DSO, especially if it is 
very compact. The material that passed during the peri-bothron 
passage the LI -point will not stay in the SgrA* Roche lobe and 
will not be accreted immediately. Only thin, hot material asso- 
ciated with a possible stellar wind of a hot and massive object 
may stay past the LI -point, but at typical escape velocities of 
such volatile material of a few 100 km/s it would have left the 
star anyhow. 

Shcherbakov & Baganoff' (2010) discuss the feeding rate of 
SgrA* as a function of radius. Their modeling suggests that 
sources X3 and X7 (Muzic et al. 2010) are still in the regime 
in which most of the inflowing mass is blown away again (see 
their Fig. 3). This may explain in part the bow- shock structure 
of X3 and X7. During its peri-bothron passage the star S2 has 
been well within the zone in which matter of its stellar wind 
could have been accreted. However, no efl'ect on the variability 
of SgrA* has been reported. The DSO peri-bothron will be at a 
larger radius than that of S2. This implies that, if the stellar wind 
from its central star is not stronger than that of S2, no enhanced 
accretion efl'ect will result from it. If the radius dependent accre- 
tion flow is clumpy or un-isotropic the DSO may still develop 
a bow- shock like X3 and X7 and a major portion of the cold 
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a) 



b) 




c) 



Fig. 18. Sketch of the relative position and motion of the Lagrange 
point LI and the DSO. 



and dusty material that was detached from the DSO during its 
peri-bothron passage will be blown away. 



5. Summary 

We have shown that most of the L'-band excess sources are in- 
deed associated with stellar objects that can be identified as H- 
or K^-band sources with the same proper motions as the L'-band 
objects. While this is the fact for many of the infrared excess 
sources we found in the inner few arcseconds of the GC, it is 
also true for all the objects that constitute the small association 
of co-moving IRS13N sources (sources JS to 7] in the nomencla- 
ture of Eckart et al. 2004). The colors of the excess sources in 
the field (excepting D3 and the DSO) have brightnesses and col- 
ors that agree well with those of the narrow line Wolf Rayet star 



candidates found in the central stellar cluster by Moultaka et al. 
(2005). 

The IRS13N sources have visual extinctions that are higher 
than those of the excess sources in the field by up to 30 mag- 
nitudes. The integral field spectroscopy shows that the IRS13N 
sources are located at a minimum of the CO-Band head absorp- 
tion in that region. This makes it unlikely that they are dusty 
late type stars older than a few Myrs. Given the colors and ex- 
tinctions of the IRS13N sources (fi to r]) they are the best can- 
didates for being young stellar objects in the overall GC stellar 
cluster with ages less than 1 Myr. They show no evidence for 
CO-bandhead absorption or emission and are then likely to be- 
long to the Class I sources of spectrally rather featureless dusty 
young stars. 

Eor the fast moving dusty S-cluster object (DSO) reported by 
Gillessen et al. (2012) we could provide a K^-band identification 
as well. While the H-K^ colors point at a significant contribu- 
tion of dust there is also the possibility that the DSO is a dust 
enshrouded stellar source and not a core-less cloud of gas and 
dust. The fact that the DSO may be a stellar source and only 
for a few months close to the SMBH has a significant influence 
on the possible disruption of its dusty envelope or disk. While 
the dusty material will certainly experience the gravitational pull 
by the SMBH SgrA* it may be rapidly inclosed again by the 
DSO's Roche lobe a few months after the peri-bothron passage 
(Eig.[T8l). While this is a possible scenario for the cold and dusty 
material any fast wind possibly associated with the DSO if it is 
a more massive, hot object, may be accreted. However, the peri- 
bothron passage distance to SgrA* is more than twice as large 
than that of S2 which did not have any reported efl'ect on the 
variability of SgrA*. It may, however, not be excluded that if 
the DSO passes by SgrA* it may develop a bow-shock similar 
to sources X3 and X7 (Muzic et al. 2010). Eurther monitoring 
of the trajectory and size evolution of the DSO as well as the 
variability of SgrA* are essential to get a complete picture of 
the scenario. We think it is currently not possible to reliable de- 
termine if and how much gas will be accreted. It is unclear if 
the object is already quite compact (as indicated by the L-band 
data) due to the influence of earlier passages. It is unclear what 
the interaction with a wind from SgrA* is and how much of the 
material will be blown away. New and detailed simulations (be- 
yond those given by Burkert et al. 2012) will be needed that 
include the DSO identification as a stellar system and the nature 
of SgrA*. 

If the IRS13N sources, in addition to the infrared excess 
sources including the DSO, are young stellar objects, this im- 
plies that the conditions of star formation at the center of the 
Milky Way may support a significant constant star formation rate 
in the overall central stellar cluster in addition to larger star for- 
mation events (in an accretion disk; Levin & Beloborodov 2003, 
Nayakshin & Sunyaev 2005, Nayakshin et al. 2006). 
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Appendix A: Appendix 

In Figs. lA.lllA.2l and lA.3l we show the relative positions of the 
stars S23, S57, S54, S63 and the dusty S-cluster object (DSO) 
over the years 2002 - 2012. in comparison with an K^-band 
image. We show the results of the model calculation (left) with 
the location of the sources on their sky projected track (dot 
interval 0.5 years; 1999-2006.5 in blue, 2007-2012 in red), an 
image of the model with the sources indicated by red circles 
(middle; image size is 380masx380mas), and a deconvolved 
K^-band image with the modeled image section overlay ed 
(right). For the years 2006/7 to 2012 a source can clearly be 
identified at the L'-band position of the DSO. The position of 
SgrA* is indicated by a red asterix. In the years 2002 to 2005 
the DSO is confused with S63. In 2006-2008 it emerges from 
the confusion and shows up between the stars S23, S63, and S54 
just north of S63. For 2009, 2011, 2012 the m^=18.9 source 
cannot be clearly identified. This is to a large part due to the 
stronger crowding at the center and to the limited data quality. 

In Fig. IA.4I we show a comparison of results from diff'er- 
ent data reduction methods and diff'erent K^- and L'-band 
datasets for 2012. In both bands diff'erent methods result in a 
clear detection of a K^- and L'-band counterpart of the DSO 
(see Ott, Eckart & Genzel 1999 for a detailed comparison of 
deconvolution algorithms in the GC field). As the DSO is close 
to the confusion limit the peak of the flux distribution within the 
diff'raction limit results in small deviations in the apparent peak 
position. 

In Fig. IA.5I we show a comparison of results from difl'er- 
ent data reduction methods and difl'erent K^-band datasets for 
2008. At both (selected as an example) epochs the difl'erent 
methods result in a clear detection of a K^-band counterpart of 
the DSO (see also Ott, Eckart & Genzel 1999). 

In Fig. IA.6I we show the L'-band identification of the DSO 
obtained via high-pass filtering. Here the source can always 
clearly be identified. For the year 2010 we have no high quality 
L'-band data at hand. 

In Fig. IA.7I we show the R.A. and Dec. values of the L'- 
and K^-band identifications of the DSO as a function of time. 
The Icr deviations between the expected positions and the 
identifications range between ±13 mas and ±25 mas (i.e. a 
quarter of the K^- and L'-band beam). 

In Fig. IA.8I we show histograms of the deviation from the 
proposed orbit (Gillessen et al. 2012) for all bands. The Icr 
deviations between the expected positions and the identifications 
range between ±13 mas and ±25 mas. 

In Figs. IA.9I we show the 1.2" angular resolution 1.3mm 
CARMA map of the mini- spiral presented by Kunneriath et al. 
(201 1) in order to outline the continuum emission from IRS13N. 

In Figs. lA.lOl lA.llI and IA.12I we present all K^-band proper 
motions of the 34 sources in and close to the IRS13E and 
IRS13N field. Off'sets in R.A. and Dec. in milliarcseconds are 
plotted against the time in years. Source labels are given as 
indicated in Fig. [2l and Tab. [3] 



field around SgrA* as shown in the overview map in Fig. [T] and 
fisted in Tab.lf The 2004 map sizes are 0.932"x0.932" (70x70 
pixels) in the K,-band and 0.945"x0.945" (35x35 pixels) in 
the L'-band. We have indicated the sources with arrows and 
labels. The grey scaling is given in bars and the contours are 
logarithmically spaced and given for display purposes only. The 
DSO is visible at the lowest contour levels in the appropriate 
plot in Fig. IA.13[ In Fig. IA.141 we marked the region in which 
D7 is located with a dashed circle. Some contour excursions 
are visible, however, now clear detection of a compact source is 
possible. For source D4/X7 see Muzic et al. (2010). 

In Figs. IA.15I and IA.16I we present all proper motions ob- 
tained for the infrared excess sources in the central 2" close to 
SgrA*. Off'sets in R.A. and Dec. in milliarcseconds are plotted 
against the time in years. Source labels and observing bands are 
given as indicated in Fig. [T] and Tab. IH The source positions 
agree to within better than about ±25 mas i.e. a quarter of the 
L'-band beam. Within the uncertainties the motions obtained at 
diff'erent wavebands agree. For those sources for which S-star 
identifications are given the positions and motions agree on 
average to within about 30 mas and 150km/s. 



In Figs. IA.13I and IA.141 we compare the K^- (left) and L'- 
band (right) identifications of L'-band excess sources in the 




Fig. A.l. K,-band identification of the DSO. 
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Fig. A.2. K^-band identification of the DSO - continued. 
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Fig. A.3. Ky-band identification of the DSO - continued. 
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Fig. A.4. Retreaving the DSO in 2012 using different methods. 
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Fig. A.5. Retreaving the DSO in 2008 using different methods. 
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201 z 



Fig. A.6. L'-band identification of the DSO. 
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U and Ks band proper motions of the DSO 
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Fig. A.7. L'- and K^-band coordinates of the identification of the DSO shown separately for R.A. and Dec. 
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Fig. A.8. Histograms of the deviations from the proposed orbit (Gillessen et al. 2012) for all bands. The deviations between the expected positions 
and the identifications in the diff'erent bands is on average less than +25 mas (i.e 1/4 of the diffraction limit in the L'-band). 
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Fig. A.9. The mini-spiral region image obtained from 1.3 mm CD array configuration with a synthesized circular beam size of 1.2''. Contour 
levels are 0.015, 0.02, 0.03, 0.04, 0.05, 0.1, 0.6, 0.9, 1.2, 1.5, and 2.5 Jy/beam. Solid black lines and a thin cyan ellipse (in the inset) indicate in the 
map and the inset the region in which the IRS13N sources are located. 
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Fig. A.IO. Proper motions of sources in the IRS13E and IRS13N field. Off'sets in R.A. and Dec. in milliarc seconds are plotted against the time in 
years. Source labels are given as indicated in Fig. [2l and Tab. [3] 
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Fig. A.ll. Fig. lAlOl continued. 
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Fig. A.12. Fig. lAlOl continued. 
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Fig. A.13. Comparison of the K^- (left) and L'-band (right) identifications L'-band excess sources in the field around SgrA*. The arrow color 
coding for the observing bands and the proper motions are given between the images in the lower part of the panel. 
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Fig. A.14. Fig. lAlSl continued. 
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Fig. A.15. Proper motions of infrared excess sources in the central 2" close to SgrA*. Offsets in R.A. and Dec. in milliarc seconds are plotted 
against the time in years. Source labels and observing bands are given as indicated in Fig.[T]and Tab.|4l 




Fig. A.16. Fig. lA.15l continued. 
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Table A.l. The L'-, K^, and H-band positions of infrared excess sources in the central 2" field close to SgrA*. Typical uncertainties 
are ±0.030". For D3 and D5 no reliable H- and K^-band identification could be given. D4 may be associated with S50 (Muzic et al. 
2008). The adopted combined relative positions are printed in bold face. 
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Table A.2. The L'-, and H-band proper motion velocities of infrared excess sources in the central 2" field close to SgrA* 
including results from a linear fit to the DSO data. The adopted combined velocities are printed in bold face. For D3, D4 and D5 no 
reliable H- and K^-band motion could be given. 
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Table A.3. H-band positions and proper motion velocities for S-source presented by Gillessen et al. (2009) calculated for the epoch 
2006.5 and assuming that 1 mas/year corresponds to 40 km/s at the distance of the Galactic Center. Following the information given 
by Gillessen et al. (2009) the typical uncertainties are ±1 mas for the positions and ±6 km/s for the velocities. For other S-sources 
used in the modeling presented in Figs. lA.l[|A.2l and rA.3l see Tab.[6j 



